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Bose Einstein Condensates, with their coherence properties, have attracted wide interest for their 
possible application to ultra precise interferometry and ultra weak force sensors. Since condensates, 
| unlike photons, are interacting, they may permit the realization of specific quantum states needed 

. as input of an interferometer to approach the Heisenberg limit, the supposed lower bound to pre- 

cision phase measurements. To this end, we study the sensitivity to external weak perturbations 
of a representative matter-wave Mach-Zehnder interferometer whose input are two Bose-Einstein 
condensates created by splitting a single condensate in two parts. The interferometric phase sensi- 
tivity depends on the specific quantum state created with the two condensates, and, therefore, on 
the time scale of the splitting process. We identify three different regimes, characterized by a phase 
sensitivity A9 scaling with the total number of condensate particles N as i) the standard quantum 
limit AO ~ l/N 1/2 , ii) the sub shot-noise AO ~ 1/N 3/4 and the iii) the Heisenberg limit A6 ~ l/N. 
However, in a realistic dynamical BEC splitting, the l/N limit requires a long adiabaticity time 
scale, which is hardly reachable experimentally. On the other hand, the sub shot-noise sensitivity 
pH ■ AO ~ 1/iV 3 / 4 can be reached in a realistic experimental setting. We also show that the 1/N 3 ^ scal- 

ing is a rigorous upper bound in the limit N — > oo, while keeping constant all different parameters 
of the bosonic Mach-Zehnder interferometer. 
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■ I. INTRODUCTION 

o 

In the last few years theoretical and experimental efforts have been devoted to the realization of an ultraprecise 
quantum interferometer (for a review see [lj h High resolution phase measurements find applications in the detection 
of ultraweak forces, as gravitational waves 0, and inertial forces The Mach-Zehnder (MZ) interferometer is a 
prototypical apparatus employing both optical and matter waves. In the optical MZ, the phase sensitivity depends 
on the intensity of the laser field, which corresponds to a limit on the average number N of photons. Typically, 
the sensitivity is bounded by the shot noise limit l/y/~N, which can be reached with classical states (for example 
, vacuum plus coherent state) as input of the interferometer. This limit, however, is not fundamental, and it can be 
surpassed with the use of nonclassical states exploiting, in this way, quantum correlations. The insuperable lower 
bound to precision phase measurements is believed to be given by the Heisenberg limit l/N Q. Its achievement 
has been theoretically demonstrated in a large body of literature with a variety of states and optimal performances 
HL IE S H> El> EH • Bose Einstein Condensates (BEC), with their coherence properties, have attracted a wide interest 
for their possible application in ultraprecise interferometry 11], and ultraweak force sensors 12. 13, 14]. Condensates, 
unlike photons, are interacting, and this property seems very promising for the realization of quantum states to use as 
input of an interferometer to approach the Heisen berg limit. For example, the recent experimental creation of Fock 
states and of very stable double- well traps [Tfl Il7| bodes well for the future of matter- wave interferometry. 

In this paper, we analyze a BEC Mach-Zehnder interferometer. The initial state configuration is prepared by 
trapping a condensate in a double- well potential with an interwell barrier large enough to create the two independent 
condensates that feed the interferometer. The height of the potential barrier is decreased instantaneously, and a 
tunneling between the two condensates is allowed for a time t_/2 = n/2Q, where f2 is the condensate tunneling rate 
between the wells. The barrier is then increased again in order to have a negligible tunneling rate. During this time, 
the interaction of a weak force with the condensates will shift their relative phase by an amount proportional to 
the energy gradient induced by the external field and the time of exposure. After a second tt/2 pulse, the relative 
number of particles is measured, and information on the phase shift is recorded. Different measurement schemes 
have been proposed, based on a positive operator value measurement |Tsj , on a p arity measurement |Tflj , on a relative 
number fluctuation measurement 0, |2(j , and on a collapse and revival detection plf . We study in detail the splitting 
processes by analyzing the produced quantum states, and giving a prediction on the sensitivity of the interferometer 
by using an error- propagation formula. We assume, for simplicity, lossless devices (beam splitters). However, losses 
of atoms can degrade the sensitivity back to the shot noise limit |22j |. With this setup we expect an improvement 
of phase sensitivity, reachable with classical states, toward the quantum Heisenberg limit. In fact, as analyzed by 
Jaaskelalnen, et al. p3|, for a complete adiabatic split, a repulsive-interaction condensate will end in a Fock state 
h/j) = \N/2)\N/2). With this state as input, the MZ phase sensitivity is expected to scale at the Heisenberg limit 
|2|. In this paper, we point out the existence of three different regimes, which, in analogy with three corresponding 
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regimes existing in the dynamical Josephson effect |29|, we call Rabi, Josephson and Fock. These three regimes are 
characterized by a MZ phase sensitivity scaling as TV" 1 / 2 , N~ 3 / 4 and TV -1 , respectively, where N is the number of 
particles in a single interferometric experiment. However, in a realistic dynamical splitting process, the 1/N limit 
requires a very long adiabatic ramping of the potential well. Its achievement is therefore strongly limited by the finite 
lifetime of the condensates. On the other hand, a sub shot-noise scaling 1/N 3 / 4 can be actually reached. We restrict 
our discussion to a two- mode analysis, using a combination of numerical and analytical tools. 
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II. TWO MODE APPROXIMATION 



In this section, we review a two-mode analysis of two weakly interacting BECs. The second quantization Hamiltonian 
of a system of dilute bosons is given by 

H(t) = Jdzi>\z)^-^^ + V(z,t)^(z) + ^Jdz¥(z)¥(z)*(z)*(z), (1) 

where ^P(z) is the bosonic field operator, and V(z, t) is the time-dependent external symmetric double-well potential, 
and g = 47r ^ - is the strength of the interparticle interaction, with a being the s-wave scattering length. The two-mode 
ansatz reads 

= Mz)a + M*)b, (2) 

where ip a ,b{ z ) can be constructed as sum and difference of the first symmetric and antisymmetric solutions of the Gross- 
Pitaevskii equation of a double well trap, correspondently to the ground and the first excited states, respectively. The 
operators and b' (a, b) create (destroy) a particle in the modes a, b, respectively. In the two-mode approximation 
the Hamiltonian of the system becomes |2J, l2fij 

H=— (a)a)aa + b r b r bb) - K(t) (a r b + Pa) . (3) 



where 

E c = 2g I dz |^ a (z)| 4 = 2g I dz \^ b {z)\\ (4) 



K(t) = - dz 



h 2 dr a {z)dMz) 



2m dz 



dz 



r a (z)v(z,t)Mz) 



(5) 



are the "one-site energy" and the "Josephson coupling energy", respectively. The operator N = h a + hb = aJa + wb 
is the total number of particles and commutes with H. By ramping the potential wells, K(t) decreases with the 
decreasing of the overlap between the wave functions. For a linear ramping, a WKB approximation [2f| gives an 



K(0)e 



-t/r 



where the effective ramping time r = 



At, 



depends on the real ramping 



exponential decrease K(t) — ii^i, , nu^n, ^^v±,^ n.m F ii, 5 . — 
time Ai ramp , on the final distance d between the wells, the height of the potential barrier Vo, and the chemical 
potential [i . Such a time-dependent configuration has been realized by the MIT group [l6j . 

We study Eq. in a phase-states representation . We write a general state in the Hilbert space of the two- mode 
system as 



+7T 



(6) 



where <fi is the relative phase between the two modes, and 



N/2 



a in<t> 



\N/2-n)\N/2 + n) 



(7) 



are un-normalized vectors of the overcomplete phase basis, written in the relative number of particles n. In this 
representation the action of any two-mode operator applied to can be represented in terms of differential operators 
acting on the associated phase amplitude ^{<f>, t). The main consequence of the overcompleteness is the non-standard 
inner product between phase vectors J7J) ( 4> \ 9 ) = j^j cos N (^j^) , which affects both the inner product between states 
© and the mean values of observable. By applying the Hamiltonian © on the states © we obtain 



H\i>) = 



(8) 



where the effective Hamiltonian is 



Heff(4>,t) 



E 3 2 K 2 (t) 



0) 
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In the first part of this work, we solve numerically and, in some limit, analytically the eigenvalue equation 

H elf {4>)^ gs {4>) = E gs ^ gs {ct>), (10) 

for different values of K . This will provide the limit of an adiabatic splitting of the condensate. In the second part 
we study the dynamical equation 

ih^^- = H ef/ (4>,tM4>,t) (n) 

for different ramping times of the interwell barrier. In both cases we give predictions of the Mach-Zehnder interfer- 
ometer phase sensitivity, as discussed in the following section. 



III. MACH-ZEHNDER INTERFEROMETER 

For a compact analysis of the Mach-Zehnder interferometer, we introduce the Hermitian operators 0, l^} 

J x = ~(tfb+&a) (12) 
4 = i(at£-Sta) (13) 
J z = \(tfa-b%). (14) 

These operators form a SU(2) Lie algebra [Ji, Jj] — itij.kJk and commute with the total number of particles N. The 
action of the Mach-Zehnder interferometer elements (beam splitter and phase shifter) on the vector J = (J x , J y , J z ) 
can be represented by rotations in three-dimensional space. In particular, the whole interferometer can be represented 
by a rotation of J around the y-axis of an angle 9 corresponding to the relative phase shift in the two arms. In our 
model, the information on the phase shift is inferred from the measure of the total N and the relative mean (J z out ) 
number of atoms in the two final condensates. We have 

(J° ut ) =-sin9 (J x ) + cos 6 (Jz) (15) 



(AJ° ut ) 2 = sm 2 9{AJ x ) 2 + cos 2 9 (AJ Z ) 2 + sinflcosfl (2(J X )(J Z ) - (J X J Z ) - (J Z J X )), (16) 

where the expectation values are taken on the input state. The sensitivity of the interferometer for the measure of 
the relative phase can be calculated by using the equation p| 

(A0) 2 = (Aj/ " t)2 . (17) 

\d(j z mt )/de\ 2 

For the optimal working pointy 9^0 and for perfect symmetric splitting of the condensate, ( J z ) — 0. By using Eqs. 
(UniEEJ), we have 

(A9f = ^0. (18) 

This equation allows us to predict the sensitivity of the MZ interferometer by knowing the input states. Because 
of its simplicity, this equation has been widely used in the literature to give predictions on the phase sensitivity of 
the Mach-Zehnder interferometer (see, for example [3lJ and ref.s therein). We note, however, that this equation is 
based on the assumption that phase distributions are Gaussian, a property that, in general, is not satisfied for states 
which are not coherent. According to the central limit theorem, to obtain a Gaussian phase distribution, we have 
to combine several independent measurements p, each having N = Nt/p particles. Therefore, equation 1)17(1 usually 
gives the right sca ling but a wrong prefactor. The exact MZ phase sensitivity can be obtained only with a rigorous 
Bayesian analysis |28j, which is rather cumbersome, and it will not be attempted here. We also remark that Ea. l(T7|) 
is obtained in the linear case, i.e. with non-interacting particles. This limit can be reached by switching off the 
interatomic interactions of the condensate after the creation of the initial state. In the general case, Eq. i|17|) gives a 
lower bound for the MZ phase sensitivity. In the phase basis, we have: 

Jx |Vj> = / + S [ sin0 ^ + (f + \) c< H |0> ' (19) 
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(20) 



and 



Jz \i>) = 



Jl M = - 



d0 <9 2 
2^ 



(21) 



(22) 



We calculate the expectation values as in equation l]18[l. taking into account the overcomplete properties of this base. 
We find 



(A9f 



where 



/-J & /-J £ W> [sin + (f + 1) cos0] 



(23) 



cos 



N 



(24) 



Analytical calculations of the quantity (|23|l can be obtained in certain interesting limits. As we will discuss later, in 
the Rabi and Josephson regimes we can consider a Gaussian phase amplitude ^>(4>,t) w ^ 4(7 ' and approximate 
the inner product as cos^ ((x — 0)/2) ~ e ^ JV / 2 (x-0) _ If er <C 7r, we can extend the integral to ±oo, obtaining 



(A9f 



N 



\2 a \l+ANa 2 )) Kl+ANa 2 ) 



2 ■ 



(25) 



In the opposite limit, when the phase uncertainty is of order of 2w, we can neglect the inner product and consider the 
phase amplitude ^((f>, t) » (1 + 7C0S</>). In this limit, 7 — + 0, we have 



{A6f 



2 + 7 2 



(26) 



(A^ + 2) 2 ' 

The whole phase sensitivity range can be exploited with the variational phase amplitude ^(4>, t) k, e 7COS ^. We obtain 



2 7 /_J r d0(cos0) Ar [cos0Ii(27Cos0) -7(sin^) 2 (I o (2 7 cos0) + I 2 (27cos0))] f_J &cj> (cos 0) N Io(27 cos < 



( f+J d</> (cos 0) w [(iV + 1) cos Ii (2 7 cos </>) - 7 (sin 0) 2 (l (2 7 cos 0) + I 2 (2 7 cos 0)) 
where l v {x) are Bessel functions of first kind and degree v. 



(27) 



IV. ADIABATIC SPLITTING 

Following the notation introduced in [2jJ in the context of quantum tunneling between BECs in a two-well system, 
we can distinguish three main regimes depending on the ratio K/Ec and the number of particles N . As we will see, 
these three regimes correspond to different input states of the MZ interferometer and, in particular, to three different 
scalings of the MZ phase sensitivity with the total number of particles. 

Rabi Regime K/E c ^> N: This corresponds to a regime in which the two wells are not completely separated, 
the potential barrier is of the order of the chemical potential, and a strong tunneling exists between the two 
condensates. In this limit, we can neglect the cos</> potential term in Eq. @, and write 



HeffW 



2 d6 2 



K 2 

E r 



■ cos 20. 



(28) 
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We note that this effective Hamiltonian does not depend on the total number of particles. A simple harmonic 
oscillator estimation gives 

which means that the phase dispersion is very small (cr? <C 1/A) and independent of N. The strong tunneling 
characterizing this regime keeps the relative phase between the two condensates well defined. However, the 
ground state of the Hamiltonian 128|) is given by a function with peaks at 4> — an d 4> — ±7r- To eliminate the 
unphysical peak at = ±7r it is necessary to take into account the correction term e^ cos ^, contained in the 
full Hamiltonian Eq.JSJ. This correction term has the same width l|29[l as the ground state of the Hamiltonian 
H e ff. The corrected phase amplitude ^(<fi) can be well approximated by a Gaussian of width given by half of 
Eq. (EHJ: 

yU) = — — *—. (30) 

\ 27T 8k) 

In the Rabi regime, taking into account that Act 2 , -C 1, from equation Q3 , we obtain 

M = 7if (31) 



which corresponds to the classical shot noise limit. We note that the phase amplitude l|3U|l does not depend 
on the total number of particles. The N dependence in A8 is merely a consequence of the inner product 
B24JI . The sensitivity scaling as ~ 1/V~N is what we expect for a MZ interferometer fed by the coherent state 
|V) - (AH)- 1 ^ +St) JV |vac). 

Josephson Regime <C <C A: In this limit, we neglect the cos(20) term in the effective Hamiltonian 10, 
obtaining 



Ec a 2 

2 



Heff (<f>) = —f-^ - KN cos 0. (32) 



As in the previous regime, we can approximate the phase amplitude with a Gaussian of width 



In this case, however, the correction term e~ cos ^ gives a negligible contribution. From equation l|25[l. we obtain 

4A\ 1/4 1 



^=UJ W (34) 
We note that this result is recovered neglecting the inner product J23J. This is a consequence of the fact that 

in the Josephson regime the two functions ip g s((f>), and tpgs{4>) — f-„ dxtp gs {x) cos w ((x — 4>)/2) have the same 
width. 

Fock Regime <C -h: The effective Hamiltonian is still given by Eq. 0, but now we have a nearly-free evolution. 
In this case the phase amplitude spreads over the whole 2tt interval and we can approximate it by 

7 cos <p 

* w =7mw,- (35) 

where 7 is a variational parameter and Io(27) is the Bessel function of the first kind of degree zero. The behavior 
of 7 as a function of the parameters E c , K, A, can be found by minimizing the total energy 



J-7T 27r 



2 



E c d 

TW~ ° S( 



r (</>), (36) 



where we have not taken into account the inner product 124|) . If we impose dE/d<p — 0, we obtain the equation 
I (2 7 )Ii(2 7 ) + (7 - ^) K(27) + Io(7)I 2 (2 7 ) - 2I 2 (2 7 )] = 0. (37) 
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In the limit 7 — > 0, we are in the Fock regime, the input state is given by \tp) = (aJ) N / 2 (w) N / 2 \vac) and the 
phase amplitude becomes flat \&((/>) ~ * l/y/2n. The Fock state is characterized by a random phase [3(|. This 
property renders the Fock state not useful in Young double-slit interferometry . In fact the interference fringes 
obtained in a single experiment are washed out by statistical averaging over many experimental runs. In this 
limit, approximating Eq. 135(1 with ^(4>) ~ l + 7COS(/>, and by neglecting terms of order o(7 2 ) in ((26(1 . we obtain 



Afl = Jf, (38) 

which corresponds to the Heisenberg limit of phase sensitivity. Correcting to this equation to higher order in 
the variational parameter 7, it would be possible to exploit Eq. ((27(1 . and express the Bessel functions 1^(7), 
v = 0,1,2 with the series expansion 

+ 00 / 7 \2fc 

= ^WWTTY (39) 

The corresponding MZ phase sensitivity is shown by the blue line of Fig.Q. 

In figure we plot the phase sensitivity with N = 1000, Ec = 0.001 msec -1 , as a function of K. We emphasize 
the three different regimes for the Mach Zchnder interferometry sensitivity characterized by a different scaling with 
the number of particles. We notice that the Josephson region becomes wider by increasing the number of particles N 
and while keeping constant all other parameters. Therefore, strictly speaking, the phase sensitivity of the BEC MZ 
interferometer studied in this paper is bounded by I/N 3 / 4 in the limit N — > 00, when all other parameters are kept 
constant. 




1(T 8 1(T 6 1(T 4 10" 10° 10 2 
K (1/msec) 



FIG. 1: Color online. Mach Zehnder phase sensitivity AO as a function of K as given by Eq.lHU. Here N = 1000 and 
E c — 0.001 msec -1 . The green lines are the analytical predictions in three regimes: i) Rabi Regime, K ^> NEc, where 
A6 = l/VW, ii) Josephson Regime, Ec/N < K « NEc, where Ad = (4K / 'E c ) 1/4 1/ 'N 3/A , and iii) Fock Regime, K < Ec/N, 
where A0 = y/2/N. The blue line is given by the variational approach with the wave function l|35|l . which reduced to Eq. I|30|l 
in the Rabi and Josephson regimes. Red points represent numerical solutions of Eq. 1101 . 



V. DIABATIC SPLITTING 



We now study the dynamical splitting of the condensate by directly solving the Schrodinger equation 111(1 , focusing 
on the dephasing process |32j . In particular, here we analyze the sensitivity of the Mach- Zehnder interferometers fed 
by states created by a finite time, diabatic splitting. As shown in [33] , initially the relative phase between the two split 
BECs has a very narrow distribution, corresponding to a BEC in a coherent state. While the height of the interwell 
barrier increases, K{t) decreases, and the phase spreads in time over the whole 2tt domain. Because of the periodic 
boundary conditions, the phase distributions eventually overlap around the region <f> ~ ±7r, developing interference 
fringes with a wavelength increasing with r. Using Eqs ((181123(1 . we study how the sensitivity of the MZ interferometer 
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changes for different ramping times r. The larger r, the more adiabatic the dynamics, and the closer it will follow the 
static model analyzed above (see Fig. QJ). If the breakdown of adiabaticity occurs in the Rabi-Josephson regime, 
the phase amplitude is sufficiently narrow so that we can study the dynamics with a variational approach using a 
Gaussian function [2^,|32. We consider |3^| 

^ )= (2 M W/* e ~^ W ' (40) 
where 5(t) and a(t) are time dependent variational parameters satisfying the differential equations 

a(t) = E c a{t)8{t) (41) 



S(t) = ^% - Ec5(t) 2 -Nk{t)e-° 2 ^ - ^e-^. (42) 

The variational wave function (|40|) is in very good agreement with the numerical calculation in the regime a ^ 1. In 
figure (0) we compare the mean-square fluctuation <j(t) obtained by the variational calculation (Eqs. I|41I42II % 1. with the 
square root of the second moment of the phase amplitude obtained by the exact numerical solution of the Schrodinger 
equation 1)11(1 . The agreement is very good until the phase amplitude touches the borders ±7r. As discussed in |32j . 
this happens when a sa 1. 



e 
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FIG. 2: Comparison between the mean-square fluctuation a, as obtained by the variational equations I|41I42|I (line) and the one 
obtained by the exact numerical solution of the Schrodinger equation lilt (points). The parameters are -/V = 1000, Ec = 0.001 
msec -1 , K(0) = 100 msec -1 , and r = 10 msec. 

In a realistic experimental setup, the two condensates initially are in the Rabi regime. While increasing the interwell 
barrier, the phase amplitude adiabatically follows the decreasing of the effective potential energy in Eq. which 
is proportional to K(t) 2 . After the breakdown of adiabaticity t a d, the phase distribution continues to spread but at 
a slower rate. When the Josephson coupling energy is sufficiently small, the dynamical evolution becomes essentially 
free, and the phase uncertainty increases as a consequence of excitations of the system. In figure © we plot the 
MZ phase sensitivity (A6*) 2 (Eq. f° r different values of r; the blue line represents the adiabatic behavior. The 

dynamics are calculated with the variational wave function l|40l) . Notice that the minimum of the phase fluctuation 
is below the point of breakdown of adiabaticity. 

We now study the MZ sensitivity by parametrizing the phase error as 

( Ad ) 2 = jp- (43) 

Here we assume that the prefactor a does not depend of the number of particles AT, as we have numerically verified 
in the limit N ^> 1. In fig. Q we plot the quantity (3 as a function of t/r = \n(K(Q)/K) obtained by calculating 
the scaling between the cases iV = 10000 and N — 1000, with fixed parameters Eq = 0.001 msec -1 , and K(0) = 100 
msec -1 . The blue line represents the adiabatic behavior, and we can clearly distinguish the three regions described 
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above. The points in fig. (J2J correspond to the MZ sensitivity at the breakdown of adiabaticity. As we can see, the 
limit 1/N can be reached only with a very long adiabatic ramping, a circumstance that is strongly limited by the 
finite lifetime of the condensate. However, the sub shot-noise limit 1/N 3 ? 4 can be achieved under currently available 
experimental conditions [h1.[17|. We remark that, by increasing N while keeping the same initial conditions, there are 
two interesting effects that take place: i) the minimum of the phase uncertainty tends to coincide with the breakdown 
of adiabaticity, and ii) the increase of the phase uncertainty after reaching the minimum is rather slow. Both these 
two effects can be understood if we approximate the evolution of the phase amplitude after t a d with a free expansion 
model (k(t) = for t > t a d) 34]. In this case 

The bigger cr 2 (t a d), the slower the dynamics of the phase amplitude. In the temporal range t a d < t < 
t a d + 2a 2 (t a d)fi/ Ec, the dynamical spreading is almost frozen. On the other hand, if we increase N by keeping the 
other initial conditions K(0), Ec and t constant, the breakdown of adiabaticity occurs at larger values of cr(t a d) 
|24|. As a consequence, by increasing N the MZ sensitivity will freeze after the breakdown of adiabaticity, and the 
minimum of phase sensitivity is maintained for a longer time. 



VI. CONCLUSIONS 



We have studied the phase sensitivity of a Mach-Zehnder interferometer fed by quantum states produced by the 
splitting a single Bose Einstein Condensate in two parts. We studied the process in the two-mode approximation, 
and we analyzed the sensitivity of the MZ interferometer by projecting into an overcomplete phase basis and by using 
an error propagation formula. We have first calculated the Mach-Zehnder phase sensitivity in the adiabatic splitting 
limit. We distinguished three different regimes (Rabi, Josephson, and Fock) characterized by a different scaling of the 
phase sensitivity with the number of particles in a single experiment. While the 1/N scaling can hardly be reached 
in realistic experiments, the limit ~ 1/N 3 / 4 , corresponding to the Josephson regime, can be achieved with current 
technology pH \ll\ , offering a considerable improvement in phase sensitivity over the shot noise 1 / \fN obtainable in 
the classical limit. 

This work was supported by the Department of Energy under the contract W-7405-ENG-36 and DOE Office of Basic 
Energy Sciences. 
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FIG. 4: Color online. Plot of the scaling parameter /3 with the number of particles (see eq. I43H ). The calculation has been 
made with N = 10000 and N = 1000, and with fixed parameters Ec = 0.001 msec' 1 , and K(0) = 100 msec -1 . The blue line 
represents the adiabatic behavior; the points correspond to the minimum of MZ phase sensitivity occurring for different values 
of t (in msec). 
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